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We show that the concept of degeneracy is the key idea for understanding the 
quantum carpet woven by a particle in the box. 



1 Introduction 

Interesting structures [1, 2] emerge in the space-time representation of the probability 
distribution for a particle in the box, as shown in Fig. 1. Three explanations of these 
quantum carpets offer themselves: Interference terms in the Wigner function [3, 4], 
degeneracy of intermode traces [5, 6] and cancelation between appropriate terms of the 
energy representation [2, 7, 8]. All of these explanations are rather involved. We there- 
fore in the present paper develop a simple argument for this surprising phenomenon. 

We identify three properties of the particle in the box as the thread of the quantum 
carpets: 

(i) The probability density involves the product of two standing waves creating con- 
tributions with the sum and the difference of the wave numbers. 

(ii) The quadratic dispersion relation connecting the energy and the momentum gives 
rise to a multi-degeneracy. 



(in 



The appropriate initial conditions enhance this degeneracy. 



The paper is organized as follows. In Sect. 2 we briefly review the important formulas 
of the problem of the particle in the box. We then in Sect. 3 give a heuristic argument 
for the quantum structures. In the Appendix A we derive a summation formula which 
allows us in Sect. 4 to cast the probability density into a form which brings out most 
clearly the canals and ridges of the quantum carpets. We conclude by summarizing the 
main results in Sect. 5. 
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Figure 1. Quantum carpet woven by a non-relativistic (left) and a slightly relativistic (right) 
particle moving in a one-dimensional box. The carpet arises from the space-time represen- 
tation of the probability density. Dark areas correspond to large probability whereas light 
areas represent low probability. We propagate a Gaussian wave packet according to the 
Schrodinger equation with the non-relativistic Hamiltonian H nr — p 2 /(2M) or with the Hamil- 
tonian H T = H n r[l — H nr /(2Mc )] approximating the relativistic Hamiltonian. In the latter 
case the straight canals and ridges of the non-relativistic box problem are curved. In both cases 
we use a Gaussian wave packet of width Ax — 0.03 L and average wave number k = 10/ L lo- 
cated at x — L/2. In the slightly relativistic example the ratio of the ground state energy of 
the non-relativistic box and the rest mass is q = 10~ 6 . 



2 The particle in the box: Fundamentals 

In the present section we summarize the essential ingredients of the problem of the 
particle in the box. In particular, we focus on the energy representation of the wave 
function. 

The probability amplitude ip(x,t) to find the particle of mass M at time t at the 
position x in the box of length L reads 

ipix^) = ^2ip m u m (x)enp (--E m t\ . (1) 
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Here the quantities 



dx <p(x) u m (x) 



(2) 



are the expansion coefficients of the initial wave packet i[)(x, t = 0) = tp(x) into the 
energy wave functions [9] 



u m (x) = 



1 f 



2L i 



— e 



ih m x 



with wave numbers 



and eigen energies 



k m = mk\ = m- 



{hk m ) 2 2 2t . 27T 

= 2M = m = m = m T ' 



In the last step we have introduced [10, 11] the revival time 

_ 4ML 2 



(3) 
(4) 

(5) 
(6) 



at which the wave function is identical to its initial form at t = 0, that is ip(x, t = T) = 
i>(x,t = 0). 

We conclude this summary of the important equations by deriving a representation 
of the probability amplitude ijj that is slightly different from Eq. (1). In Sect. 4 we will 
use this expression to bring out the location and shape of the structures in the quantum 
carpets. 

We substitute the energy wave function u m into Eq. (1) and use the expressions 
Eqs. (4) and (5) for the wave numbers k m and energies E m . We then arrive at 



1 °° f 



x 2t 
i-nm | — — m— 
L T 



- tp m exp 



x 2t 
- ' rn n \ — + m— 



When we define the expansion coefficients ip m for negative values of m by 

1p-\m\ = ~1p\m\ 

we find the compact representation 



v m— — oo 



r. fx 


2A" 


mm 1 — 









(7) 

(8) 

(9) 



of the wave function. 
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3 Quantum carpets: A heuristic argument 



In the present section we use the expression Eq. (1) to show that the probability density 
consists of four terms: Two terms correspond to the classical trajectories, whereas the 
other two represent the striking canals and ridges of Fig. 1. 
Since the structures appear in the probability density 

W{x,t) = ^*{x,t)i/){x,t) (10) 
we now use the energy representation Eq. (1) of to find W and arrive at 

M 



W(x,t) 



l 7f E C^f^ 1 - e~ tk ™ x )(e tk " x - e~ ik " x ) exp 



2L 



m,n—l 



2M 



(11) 

Here we have used the expressions, Eqs. (3) and (5), for the energy wave functions and 
the energies. 

When we multiply out the individual waves in the product of the two energy wave 
functions we recognize that the probability 



r(+) 



(-) 



consists of four contributions. The terms 

1 °° { 
! q t ] ( x , t) = -— E exp | ±i(k m + k n ) 



m,n— 1 



x zh (k m k n ) 



u 

2M 



(12) 



(13) 



arise from the multiplication of the two co-propagating waves exp(±ik m x) and exp(±ik n x) 
in the two energy wave functions. Note that the relation 



(14) 



has allowed us to factor out the sum k m + k n of the wave numbers in the expression 
Eq. (13). This creates the difference k m — k n of the wave numbers in the expression in 
the square brackets. 

In contrast the terms 



r(±) 



x,t) 



1 °° { 
2L E exp I ±i ( fc ™ ~ 

m.n—l ^ 



X i (k m + fc n ) 



fit 

2M 



(15) 



are a consequence of the multiplication of the two counter-propagating waves exp(±ik m x) 
and exp(=pifc n x) in the two energy wave functions. Here the factorization property 
Eq. (14) has led to the difference k m — k n of the wave numbers outside the square 
brackets and the sum k m + k n inside. 
The phases 



4> { m}n(x, t) = j±(k m - k n ) 



— ± (m — n)-J— 
L 'T/2 



(16) 
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and 



t) = j±(k m + k n ) JjL. = ^±( m + n)^ (17) 



in the square brackets of Eqs. (13) and (15) correspond to straight lines in space-time. 
The steepness of these world lines is determined by the difference and the sum of the 
quantum numbers m and n. As it was shown in Ref. [5], this opens the possibility for 
a multi-degeneracy: Different pairs of quantum numbers m and n can give rise to the 
same difference m-nor sum m + n. Therefore many world lines can lie on top of each 
other enhancing in this way the contrast of the structures. 

It is the expansion coefficients ip m that decide the question of enhancement or sup- 
pression of these world lines. To understand this we consider a distribution ip m of energy 
excitations which satisfies the factorization property 



Hence we can replace the product V'mV'n 01 the initial expansion coefficients by another 
product of new expansions coefficients and ipi ■* which now only depend on the 
sum s and the difference r of the quantum numbers. Any Gaussian wave packet satisfies 
this condition. 

We conclude this section by considering a Gaussian wave packet centered at quantum 
number m and width Am such that 1 <C Am <C m. In this case we find a clear 
separation of classical and quantum trajectories contributing to the probability density 
[5]: The terms 1^ contain the classical trajectories whereas the terms Iqt^ are the 
origin of the carpet. In order to bring this out we recall that the terms 1^ contain 

the phase <&m}n with the sum of the quantum numbers. Due to the Gaussian weight 
factor with its maximum at m > 1 this gives rise to large quantum numbers of the 
order of 2m. The steepness of the corresponding world lines is therefore proportional to 
TrT x <C 1. Consequently the world lines are rather flat and correspond to the classical 
trajectories. In contrast the terms Iqt^ contain the phase 4>m}n with the difference of 
the quantum numbers. This corresponds to steep world lines — the striking canals and 
ridges. 



4 A new representation of the probability density 

In the preceding section we have shown that for the evaluation of Iqt^ and 1^ it is 
natural to introduce the new summation indices m±n. In the present section we pursue 
this idea. However, we do not start from Eqs. (13) and (15) but from Eq. (9). 
The probability W then reads 



W{x, t) = — ex P \ -in(m - n) 



(m + n) 



(19) 
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With the help of the summation formula 

oo -. oo oo poo r -i i 

E = 2 E E (- 1 ) J 7 ^/ 2 c?+/9) '2 ( ^ p) 

m.n— — oo i= — oo j — — oo oo L 

derived in the Appendix A the probability density takes the form 



i,j=-oo 

Here we have introduced the Wigner function 



2> Xj.iO*,*) 



exp(iTrZp) (20) 



(21) 



J —CO 



cxp(-iTrp^), 



(22) 



of the expansion coefheients ^ TO . Note that ^>[/x] is a continuous extension of tp m such 
that ip\p] = f° r A* = m - 
Moreover, we have defined 



(23) 



Expression Eq. (21) for the probability density is the main result of the present 
paper. It brings out most clearly that the probability density \ip(x, t)\ 2 consists of 
a superposition of structures tyW) aligned along straight lines defined by Xj,l( x ,t). 
The Wigner function f tne 

expansion coefheients determines the shape of these 

structures. 

We conclude this section by making contact with the expression, 



w{x,t) = — e (-i) j ''* (w W*>*),Pi], 

derived in Ref. [3]. Here = jnU/(2L) and 



(24) 



$W(a;,p) = -L y°° di/^(a;+|) 0* (z - |) exp(-ipy/ft) (25) 



is the Wigner function [12] of the superposition wave function 

<p(x) = ip(x) - ip(-x) 



(26) 



built out of the original wave packet and its mirror image. 

The similarity of the two expressions Eqs. (21) and (24) for the probability density 
already suggests that the Wigner function ty( w > in (fi, £) space is related to the Wigner 
function $W in (x,p) phase space. Indeed in Appendix B we derive the relation 



-fi 



(27) 



Hence the two arguments fi and £ appropriately scaled now play the role of momentum 
p = hk^ and position x = L£. 
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5 Conclusions 



In the energy representation the probability of finding the particle at time t at position x 
is a double sum over all energy quantum numbers. The quadratic dispersion relation of 
the free particle allows us to express this double sum by another double sum extending 
over the sums and differences of the quantum numbers. In this way we represent the 
probability distribution in space-time as a superposition of structures along straight 
world lines. Their steepness and their starting points are determined by integers. The 
shape of the initial wave packet, that is the form ip m of the initial excitation of the 
energy eigenstates governs the shape of these structures. 

We emphasize that the expression Eq. (21) for the probability density derived in 
the present paper is very different from the one which follows from the representation 
[2, 10, 11] 

i/> (x,t = ±T + &t) = jr wl r) vU- l -2L,At) - £ wf ] J-x+ l -2L,At\ , 

I — — OO I — — OO 

(28) 

of the wave function ip in the neighborhood of a fraction q/r of the revival time T. Here 



(r) 



1 r_1 

;E ex p 



p=0 



-2ni (p 2 - — p- 
\ r r , 



(29) 



denotes the Gauss sums [13, 14] and 



ip(x,t)= / dx'G free (x,t\x',0)<p(x') (30) 
Jo 

is the initial wave function ip(x) = ip(x,t = 0) propagated freely according to the 
Green's function Gf ree of the free particle. 

Indeed the above formula Eq. (28) is a local representation in space-time whereas 
Eq. (21) is a global one: It depicts the probability density Wasa superposition of struc- 
tures along vertical and tilted world lines, whereas the revival representation Eq. (28) 
uses a superposition of structures along lines of constant time, that is along horizontal 
world lines. 

We have obtained these results using the summation formula derived in the Ap- 
pendix A. This formula has a much wider range of application. For example, it im- 
mediately provides an expression for the slightly relativistic particle. In this case the 
structures are not along straight but curved lines as shown in Fig. 1 (right). Moreover, 
it provides insight into the development of fractal canals discovered in Rcf. [2], when 
the initial wave packet is uniform. However, space does not allow us to go deeper into 
these topics of future publications. 
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Appendix A: A useful summation formula 

In this appendix we derive two different but equivalent representations of the double 
sum 

oo 

1= ]T f m ,n (31) 



m,n— — oo 



with coefficients f m .n- 

Our derivation relies on the introduction of the new summation indices m + n = s 
and m — n = r. Note however, that this definition puts certain restrictions on s and r. 
Indeed, we have to distinguish two cases: (i) when m and n are both even or odd, and 
(ii) when one of them is odd and the other is even. In the case (i) we find that m — n 
and m + n are both even. Hence we have the substitutions 

m — n = 2r and m + n = 2s. (32) 

In the case (ii) we find that m — n and m + n are both odd, which leads to the definition 

m — n = 2r + 1 and m + n = 2s + 1. (33) 

We therefore find the rule 

oo oo oo 

^ ^ fm.n ^ ^ fs+r,s — r "i" ^ ^ fs+r+l,s—r (34) 

m,n— — oo r,s— — oo r,s— — oo 

for replacing the original sums extending over m and n by new sums extending over r 
and s. 

We can combine the two terms in Eq. (34) into one, when we replace either the 
summation over r, or the one over s, by an integration. The Poisson summation formula 
[15] 

oo oc poo 

^2 9m= ^2 / dfj,g[fi\ exp(2nilfi) (35) 

m— — oo l= — oc ~ 00 

allows us to do this in an exact way. Here g[fi] is an extension of the function g m to the 
whole real axis such that takes on the values g m at integer values fi — m. 
When we apply the Poisson formula to the summation over r we arrive at 

OC f OO QQ 

I = \ / dpf[s + p,s-p] exp(2nilp) 

+ J2 dpf[s + p + l,s-p]exp(2mlp)\, (36) 
which after the substitutions p = 2p and p = 2p + 1 in the two integrals takes the form 



1 OC ( OO „oo 

= E E #/ 



-oo ^s— — oo 



i(2 S + / 5),i(2 S -p) 



exp(inlp) 
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+ 



E (-!)'/ 



dp/ 



I(2* + l + p),l(2« + l-p) 



exp(iwlp) } . (37) 



Here we have used for the last integral the relation 

exp(nrZ) = (-1)', 



(38) 



and have written the arguments of the function / in a way that brings out most clearly 
that the two integrals are the even and odd terms of a single sum. The only obstacle 
left before we can combine these two terms is the term ( — 1)'. When we recall that 



{-!)* 



(-l) 2sl = 1 for j = 2s 

(_l)(2«+i)i = (_i)i forj = 2s + l 



we find indeed 



^ 1 °° oo poo 

E E E(-W 

n.= — no r— — r^ i = -^ J—oo 



m,n— — oo 



i=— oo j=— oo 



(39) 



exp(z7r?p). (40) 



We conclude this appendix by presenting a different expression for the double sum / 
which follows from Eq. (34) when we replace the summation over s by an integration 
using the Poisson summation formula, Eq. (35). In this case we find following the same 
train of thought 



E /.»,» = 5 E EM) 



oo oo 



J I 



rn.n— — oo 



(= — OO J — — QO 



daf 



exp(i7rZcr). (41) 



We recognize that the two representations are different: In the one of Eq. (40) the 
integration variable p enters in an asymmetric way whereas in the one of Eq. (41) the 
integration variable <r appears in a symmetric way. Nevertheless, both representations 
are completely equivalent. 



Appendix B: Relation between the two Wigner func- 
tions 



In this appendix we relate the Wigner function 



cxp(— inpt;) 



of the continuous extension ip[p] of the expansion coefficients 



= / dx(p(x) u m (x) 
Jo 



(42) 



(43) 
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to the Wigner function 

<& w (^)^/_°°^(* + f) **(*-!) e ~ w?l (44) 

in position a; and momentum p of the superposition state 

</>(#) = <p(x) — ip(—x). (45) 

For this purpose we first note that the initial wave packet vanishes at the walls at 
x = and x = L. Moreover, it vanishes outside of the box. We can therefore extend 
the integral in the definition of the expansion coefficients to — oo and +oo, that is 



2 '■ ' 
L 



dxip(x) sin (mir^j 



(46) 



where we have used the definition Eq. (3) of the energy eigenfunctions. 

When we substitute this expression into the definition of the Wigner function Eq. (42) 
we arrive at 



1 /■ OO /• oo 

* W (M,0 = -2lJ dx 'J dx"<p'&)<p(x") 



x < exp 



x' + x" 



exp nr/i 



— exp y—mfi 
+ exp I — iir[i 



L 

x' + x" 



dp exp 

' — oo 

/CO 
dp exp 
-co 
poo 

dp exp 
dp exp 



mp 



2L 



tirp 



2L 



iirp 



iirp 



2L 
x'-x' 



2L 



(47) 



We now introduce the new integration variables x" = —x" in the second, x' = —x' 
in the third and the pair x' = —x' and x" = —x" in the forth term of the brackets. 
This allows us to combine the four contributions which yields 



dx 



dx" [<p*(x') - ip*(~-x r )} [<p{x") - <p(-x")] 



X CXp ITTfl 



2L 



a- 



Here we have also made use of the relation 

dpe™ pc = 2 5(Q. 



I 

J — ! 



(48) 



(49) 
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When we now perform the integration with the help of the delta function we find 

^> {W) (jt, = -2 J °° dx" 4>* (x" + 2L£) 4>(x") exp ^2tt^ — ± ^ (50) 

which with the new integration variable —y/2 = x" + L^ reads 

* W (M, = -J°° dy 4>* fa - |) - |) exp(-iir W /L) (51) 



or 



*W(^,0 = /°° + exp(-i7r W /L) (52) 

where we have made use of the symmetry relation (j){—x) = —(f>(x). 

When we compare this expression to the Wigner function & w ">{x,p), Eq. (44), in 
position x and in momentum p we find the relation 

*WM = toh*w(Lt,^. (53) 

We conclude this appendix by noting that the operation of complex conjugation 
enters differently in the two Wigner functions: In (x,p) phase space we take the complex 
conjugate of the wave function with the argument x — y/2. In contrast, in the Wigner 
function in (jj,, £) space we take the complex conjugate of the expansion coefficient with 
the argument u, + p/2. 
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